Quantum computations with atoms in optical lattices: 
marker qubits and molecular interactions 
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We develop a scheme for quantum computation with neutral atoms, based on the concept of 
"marker" atoms, i.e., auxiliary atoms that can be efficiently transported in state-independent pe- 
riodic external traps to operate quantum gates between physically distant qubits. This allows for 
relaxing a number of experimental constraints for quantum computation with neutral atoms in mi- 
croscopic potential, including single-atom laser addressability. We discuss the advantages of this 
approach in a concrete physical scenario involving molecular interactions. 
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I. INTRODUCTION 

Manipulation of cold atoms in microscopic traps is one 
of the major highlights of the extraordinary progress 
experienced by atomic, molecular and optical (AMO) 
physics over the past few years, and has lead to impor- 
tant successes in the implementation of quantum infor- 
mation processing 0. By employing a quantum phase 
transition it is possible to load large numbers of neutral 
atoms in highly regular patterns within an optical lattice 
. This system is very promising both in terms of quan- 
tum simulation of condensed matter physics, and more in 
general of quantum information processing. Hence, over 
the last few years, several implementations of neutral- 
atom quantum computing, exploiting various trapping 
methods and entangling interactions, have been propose d 
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In this paper we study quantum computing with neu- 
tral atoms in optical lattices based on the concept of 
marker and messenger atoms. We consider a situation 
where qubits are represented by the internal longlived 
atomic states, and these qubit atoms are stored in a 
(large) regular array of microtraps realized by an opti- 
cal lattice. These qubit atoms remain frozen at their 
positions during the quantum computation. In addition 
to the atoms representing the qubits, we consider an aux- 
iliary "marker atom" (or a set of marker atoms) which 
can be moved between the different lattice sites contain- 
ing the qubits. The marker atoms can either be of a 
different atomic species or of the same type as the qubit 
atoms, but possibly employing different internal states. 
These movable atoms serve two purposes. First, they al- 
low addressing of atomic qubits by "marking" a single 
lattice site due to the marker-atomic qubit interactions: 
this molecular complex can be manipulated with a laser 
without the requirement of focusing on a particular site. 
Second, the movable atoms play the role of "messenger" 
qubits which allow to transport quantum information be- 
tween different sites in the optical lattice, and thus to 



entangle distant atomic qubits. 

The first key element in our scheme is the transport of 
marker (or messenger) atoms in an off-resonant time- 
dependent superlattice. By changing laser parameters 
with an appropriate protocol we move the marker atoms 
from site to site while leaving the qubit atoms frozen 
at their respective positions. We note that to move a 
marker atom only the global laser parameters generat- 
ing the superlattice need to be changed. In the case of 
several marker atoms on a lattice (arranged e.g. in a 
certain spatial pattern) they will be moved in parallel 
by these global lattice operations. The time scale for 
these lattice movements can be of the order of the os- 
cillation period in the confining lattice potential. In ad- 
dition, two distinctive properties of the scheme are: (i) 
The superlattice can be realized by a very far-offresonant 
optical lattice. Thus there is no requirement for a qubit- 
(or spin-dependent) optical lattice as in the case of colli- 
sional gates which in case of Alkali atoms require tuning 
of the lattice laser between the excited atomic fine struc- 
ture states. This allows to strongly suppress decoher- 
ence due to spontaneous emission in the present scheme, 
(ii) There is significant freedom in choosing the internal 
atomic states representing the qubits: in particular, we 
can choose atomic states corresponding a "clock transi- 
tion". These clock states are insensitive to the (stray-) 
magnetic fields, again improving decoherence properties 
of the atomic qubits. Also this is in contrast to mov- 
ing atoms in spin-dependent lattices for collisional gates, 
where the qubit states are typically very sensitive to mag- 
netic fields. 

A second key element is that we employ resonant 
molecular interactions between marker and qubit-atoms, 
as provided by magnetic or optical Feshbach resonances. 
This implies two features of the present scheme: (i) Due 
to the resonant character combined with the spatial con- 
finement of atoms in the optical lattice, these interactions 
can be comparable to the trap spacing in the optical lat- 
tice, and thus the time scale of operations becomes of the 
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same order of magnitude as the one for the transport in 
the lattice, (ii) In addition, these resonant molecular in- 
teractions can be made internal state (qubit) dependent 
which gives a mechanism for entangling the marker and 
atomic qubits, and to perform swap operations of the 
atomic qubit to the marker atom. 

The article is organized as follows: in Sec. |H] we in- 
troduce the general concept of quantum computing via 
"marker" qubits, and we specialize it to the case of 
atomic qubits in optical lattices. In Sec. 11111 we develop 
and simulate a procedure to effect selective atom trans- 
port in spin-independent lattices. Sec. IIVI describes the 
theory of resonant collisions in confined geometries, suit- 
able for the treatment of Feshbach resonance in tightly 
confining traps. Sec. E] discusses the dynamics of one- 
and two-qubit operations using the aforementioned in- 
gredients. Conclusions are drawn in Sec. IVII 



II. CONCEPTS OF QUANTUM COMPUTING 
WITH "MARKER" ATOMS 

A. General concept 

The scheme we are introducing is based on a quan- 
tum register formed by separately stored qubits, that 
never interact directly with each other. To mediate en- 
tangling operations between different register (i?-type) 
qubits, we introduce "marker" or "messenger" (M-type) 
qubits, that can be transported through between differ- 
ent register locations. Direct coupling can only take place 
between a register and a marker qubit. In the simplest 
situation there is only one M qubit present in a certain 
register location Ri. Different operations are then possi- 
ble (see Fig. 0: 

(i) The M qubit can be transported forward and back- 
ward throughout the string of register qubits thus 
being able to reach an arbitrary location Ri . 

(ii) A local interaction between the M and the R qubit 
may be activated to perform single and two qubit 
gates. 
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FIG. 1: Basic operations with a marker atom on a quantum 
register: (i) forward and backward transport steps, (ii) local 
interaction with register qubits. 



The role of the M qubits in our scheme is twofold. On one 
hand they will allow us to address single register atoms 



without the need for addressing single lattice sites, i.e. 
they act as a "marker" for a certain register atom. On 
the other hand they act as information carrier performing 
effective entangling operations between physically distant 
register qubits. In this case they act as a "messenger" 
transporting quantum information. However, to simplify 
language, in the following we will denote the M qubits 
always as marker qubits or marker atoms. 

The interaction we apply depends on the logical state 
of both of the involved qubits, i.e. it enables us to per- 
form two-qubit gates like a controlled-phase or a swap 
gate between the marker and the register atom. This 
in turn, combined with forward and backward transport 
operations (i), allows to construct sequences of opera- 
tions that enact two-qubit gates between arbitrary regis- 
ter qubits. This works as follows (see Fig. [2J: 

a) the state of qubit Ri is first swapped onto the 
marker M; 

b) the marker is transported to location j, where a 
two-qubit gate is performed between M and Rj] 

c) the marker is transported back to location i and its 
state is swapped back onto Ri. 



a) 



b) 



c) 
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FIG. 2: Realizing an entangling operation between distant 
atoms based on the elementary steps described in Fig. a) 
swapping the first qubit onto the marker atom; b) transport- 
ing the marker atom unto the second qubit and local interac- 
tion; c) transport back to the first qubit and inverse swap. 



At the end of the process, the marker qubit recovers its 
initial state, while the net effect is that a gate opera- 
tion has been performed between register qubits Ri and 
Rj. Beside logical operations, this allows for the creation 
of distant entangled (e.g., EPR) pairs that can be subse- 
quently used for teleportation between different quantum 
memory locations, for state purification in error correct- 
ing protocols and for scalable probabilistic gates pfl |. 
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B. Implementation of the concept with atoms in 
periodic trapping potentials 

In the following we want to briefly outline how the 
above general concept can be implemented. A detailed 
description can be found in Sec. IIIIIVI As described 
above the two key ingredients of our scheme are (i) the 
transport of marker atoms and (ii) the application of a 
strong local interaction. We concentrate in this work on 
an implementation with neutral atoms stored in a two- 
component optical superlattice. However, our scheme 
may be transferred to other systems, including atom 
chips |2l| . We consider single atoms stored in the ground 
state of separate wells which we model as a ID periodic 
potential (see Sec. IIII At . with a simple filling pattern 
of one register atom Ri every second lattice site. The 
ground states of the remaining sites may (or may not) 
be occupied by marker atoms, which can be of the same 
species as the register atoms, and the tunnel coupling 
between neighboring sites is assumed to be negligible, so 
that marker atoms do not interact with register atoms 
unless the potential is modified. The quantum informa- 
tion is stored in two appropriate internal states |0), |1) 
of the atoms (which will be specified in Sec. HV(I . 

As described in Section lHTl the transport of the marker 
atoms (i) is realized by globally changing the external 
lattice control parameters which allows for creating a pe- 
riodic array of double-well structures with different well 
depths. In this way the marker atom can be transferred 
from its initial site into the first excited state of one of the 
neighboring wells, while the register atom located there 
(as well as any other register atom in the lattice) remains 
in its trap ground state. From here the marker can be 
transported further to the ground state of the next site 
(which is not occupied by an atom) or back again to its 
initial position. By repeating these transport steps the 
marker atoms can be transported to an arbitary lattice 
site. This scheme avoids there being two permanently 
interacting atoms at any lattice site. 

Due to the fact that the lattice parameters are changed 
globally, all marker atoms undergo the same, parallel 
movement. Thus, when more ground-state marker atoms 
are introduced at different sites, a certain lattice trans- 
formation will transport all of them in the same way. By 
suitably choosing the pattern of marker atoms, multi- 
qubit operations can be carried out in parallel or with 
pre-defined patterns (an example is encoding and syn- 
drome extraction for error correction). 

When a marker and a register atom are at the same 
site we realize the coupling (ii) of Fig. ^ by making use 
of the strong molecular interaction between the marker 
and the register atom, which can be controlled by an 
external magnetic field giving rise to a Feshbach reso- 
nance (23. The physics behind this mechanism, as well 
as the gate operations, will be detailed in Sees. HVl and 
fVl respectively. Of course, this sort of interaction can 
be employed in any neutral-atom quantum computation 
proposal. In this paper, we will outline its general fea- 



tures, and we will focus on its specific use in the context 
of marker-atom quantum computing. 

The principle of the single qubit gate is the following. 
The marker atom is transported to the register atom 
we want to address. Then the molecular interaction is 
"switched on" via an external magnetic or optical field, 
i.e. we perform a Feshbach ramp which leads to a level 
splitting of the atomic states. Clearly this splitting is 
only present at the site with two atoms. With appropri- 
ately detuned external lasers we can then perform arbi- 
trary single qubit rotations 

|0) — >cos(a)|l) -isin(a)e iv |0) 

|1) — ►-isia(a)|l) + cos(a)e i¥> |0) (1) 

where the angle a is given by the interaction time with 
the lasers (and their intensities) and the phase <ft is deter- 
mined by the dynamics of the Feshbach ramp. The lasers 
have not to be focused down to the lattice constant. The 
spatial width is merely limited by the distance to the next 
marker atom (except if we want to perform the same ro- 
tation there). 

The principle of two qubit gates is already shown in 
Fig.[2]where in step b) we either perform a swap operation 
between register and marker atom, \e1e2} — > |e2£i), or a 
phase gate, \eie 2 ) — > exp[i</?(l — ei)(l — ea)]|ei£2) where 
£1.2 € {0,1}. As we will describe in Sec. [3 the phase 
gate as well as the swap gate are again based on the 
tunable molecular interaction between two atoms at one 
lattice site. In the first case the phase is acquired by 
a Feshbach ramp which affects only the state |00) while 
in the case of the swap gate we need, similar as in the 
case of the single qubit rotation, an additional laser field 
which couples resonantly the state which is shifted by the 
molecular interaction and the states 1 01) and 1 10) . The 
laser will again affect only the sites where two atoms are 
present thus again it does not have to be focused. 

Beside relaxing addressability constraints, our scheme 
bears several other advantages: for instance, it does not 
require a state-dependent lattice [2^. As described in 
Sec.|Jthat method has a couple of disadvantages and, fur- 
thermore, the realization and stabilization of such poten- 
tials poses a major experimental challenge. In our scheme 
the quantum register logical state never gets entangled 
with the atomic motion, eliminating a major source of 
decoherence. Even collisional phases, acquired by the 
marker atom while being transported over occupied lat- 
tice sites, can be made state- insensitive by an appropriate 
choice of the atomic hyperfine states (a typical example 
being Rb, for which the singlet and the triplet scattering 
lengths coincide), thus contributing only a global phase 
to the evolution of the whole register. 

III. ATOM TRANSPORT IN TIME 
DEPENDENT SUPERLATTICES 

The transport scheme we describe in this section makes 
use of a time dependent optical superlattice configuration 
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Laser pair b 



the atomic dipole moment d and the electric fields 




Laser pair a 

FIG. 3: Laser setup. The plane containing the wavevec- 
tors fci,2 of a laser pair a intersects the plane containing the 
wavevectors £3,4 of a laser pair b at the x-axis. The planes 
can have a finite angle between them. The angle between the 
lasers in each of the pairs are 9 a and Ob, respectively (see also 

H). 



which can be far-offresonant from the relevant optical 
transition to avoid spontaneous emission and is not spe- 
cialized to specific atomic species. The atom transport is 
independent of the considered internal states and allows 
for using the m = states of different hyperfine man- 
ifolds, i.e. states of a "clock transition" which are not 
affected by external magnetic fields. 

In the following we will describe the laser configuration 
which is necessary to realize the superlattice and detail 
the transport of single atoms in the periodic potential by 
changing the intensities and phases of the lasers. We will 
furthermore discuss optimization methods. 



A. Realization of the superlattice 



For the realization of the superlattice potential we pro- 
pose using a configuration of four intersecting lasers like 
it was used in [24|. The setup is shown in Fig. EI Two 
pairs of laser beams intersect with an angle a and 9 b , 
respectively The lasers of frequency to a (pair a) and u> b 
(pair b) interact with atoms which are considered as two 
level systems with transition frequency luq. In an inter- 
action picture, the Hamiltonian for the internal degrees 
of freedom of an atom can then be written as 



H = H + HA a (T + cr- 



1 



<j+d(E a {r) +^ 6 (r)e i 



(2) 



where Ho = p 2 /2m is the kinetic energy operator. Fur- 
thermore we introduced A a ^ b — luq — Lo a .b, S — uj a — uj b , 



e,e 



i(kjr+(j>j) 



Ea(r) = 

3=1,2 

E b (r) = £^> i( " V+03 \ 



(3) 
(4) 



where £*, are the normalized polarization vectors, £j the 
field amplitudes and the wavevectors have the magnitude 
k a = \h\ = \k 2 \ = oj a /c, k b = \k 3 \ = |ft 4 | = w b /c. 
The laser pairs are assumed to be far off detuned from 
atomic resonance and from each other so we can adia- 
batically eliminate the upper atomic level and obtain an 
effective Hamiltonian in position representation 



^ ^ 2 \dE a (r)f \dE b (r)\* \ 

2m \ 4hA a 4hA b ' 1 ' 



According to the geometry of the laser setup the second 
term of Eq. (JSJ) can be written as 



U(x) = - 



COS (2^X8^1(6^/2) + 01 - (f>2 ) 



4hA a ^ 1 Jl 3 4hA b 

3=1,2 b j=3,4 

(dh)(d*r 2 )£i£ 2 
4hA a 



4hA b 



cos (2k b xsm(9 b /2) + (f> 3 — (f> 4 



(6) 

= U + Ui cos (2k a x sm(e a /2) + fa - fa) 
+ U 2 cos (2k b x sin(0 6 /2) + 3 - fa) , (7) 



where Uo is merely a constant. In Sec. lIII Bl we will choose 
it according to Eq. ((Tn|) . If we set fa —fa, <p = fa ~ fa, 
9 a = ir (counterpropagating lasers) and 



2 arcsin 



k a 
2h 



2 arcsin ( — 



7T 

3' 



(8) 



the potential is of the type 

U(x) = U a + Ui cos(2kx) + U 2 cos{nx - fa. (9) 

The potential © leads to a particle confinement along 
the x axis. We will assume in the following the confine- 
ment in the transverse directions y and z to be much 
stronger than along the x-direction so that we have ef- 
fectively a one dimensional system. 

B. Single atom transport 

The transport of an atom through the lattice is 
achieved by varying the amplitudes Ui(t) and the relative 
phase fat) of the two lattice components, which is done 
by changing the intensities and phases of the lasers-see 




FIG. 4: Potential U(x,t) at two different times and V = 
100 E r (with recoil energy E r = h 2 K 2 /2m). At t = 
iti(0) = 112(0) = (dashed line) while at ti > we set 
ui(ti) = 0.5, u 2 {t 2 ) = 0.2 and a = 1, I = (solid line). 
In the latter case the difference between the minima is ap- 
proximately 20 E r and the height of every second barrier is 
reduced by approximately 50%. 



Eq. ©. The potential © becomes then time dependent: 
U(x) -► U{x,t). 

For the description of the transport process it is useful 
to rewrite the potential @ depending on two parameters 
ui(t) and U2{t): 



U Q (t) 
Ui(t) 

u 2 (t) 

0(t) 



j[2-ui(t) + « 2 (*)], 
V 



[2-tti(t)-u 2 (t)], 



2-V u i(*)+«2(t) 2 , 
(j arctan [u 2 (i) /ui (i) + Z7r] 



(10) 

(11) 

(12) 
(13) 



At time t = we set U 2 (t) = 0, i.e. ui(0) = u 2 (0) = 
and thus we simply have a cosine-potential of depth 
V and periodicity do = 7t/k. However, in general the 
periodicity of the lattice is a\ = 2tt/k and the shape 
of the optical potential can be designed in the following 
way: 

The parameter Ui(t) controls approximately the height 
V[l — ui(t)] of every second barrier depending on the 
parameter I £ {0, 1}: If I = the height of every "odd" 
barrier is changed while in case of I = 1 the height of 
every "even" barrier is changed, i.e. in the first case 
barriers with maxima at x = (2j + l)ir/K and in the latter 
case barriers with maxima at x — 2jn/K are modified 
for 1/2 = and j £ Z. By changing ui(t) we can thus 
create a specific periodic array of double well potentials. 
The parameter 1x2 (0 controls additionally the difference 
of the minima Vu2(t) of such a double well potential, 
while a £ {±1} determines if the left (cr = 1) or right 
(cr = —1) well is raised. 

An example is shown in Fig. 0] for I = 0, a = 1 and two 
different values for u\ and u 2 . The potential U(x, t) is 
given in units of the recoil energy E r = h 2 K 2 /2m where 
m is the atomic mass appearing in the time dependent 
Schrodinger equation 



2m dx 2 



(2j+1/2)jt/K x (2j+3/2)it/K 
b) 100 




(2j+1/2)it/K x (2j+3/2)jt/K 
c) 100 



e) 100 




(2j+1/2)jt/K 



(2j+3/2)it/K 



ip(x,t) + U(x,t)i/)(x,t) (14) 



FIG. 5: Illustration of steps (i)-(iv) of the transport process 
as described in the text. Shown are also the square of the 
absolute value of the single particle wavefunctions of the atom 
which is transported (i.e. the "marker atom" ) , |f/> M (2;)| 2 (solid 
lines), and of the register atom which is supposed to remain 
at its lattice site, \ip R (x)\ 2 (dashed lines). The example shown 
requires a — 1 and / = 0. 



which has to be solved for the study of single atom trans- 
port. The elementary steps of the atom transport are 
done by tunnelling in the double well potentials. An ex- 
ample is shown in Fig.[S]for a = 1, 1 = and V = 100 E r : 

At the initial time t = 0, where wi(0) = 1*2(0) = 0, 
we consider two neighboring wells, with one atom in the 
motional ground state of the left well (Fig. |SJi). The 
probability densities \4> M {x : t)\ 2 of the atom are indi- 
cated by the solid lines in this figure. The superscript 
M indicates the wavefunction of the atom to be trans- 
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ported, i.e. the wavefunction of the "marker atom" as 
introduced in Sec [n] Also shown in this figure by the 
dashed lines are the probability densities of the "register 
atom" \ij) R (x,t)\ 2 , initially located in the ground state of 
the right well in this example and which is supposed to 
remain at its lattice site during the transport. Our goal 
in the process described here is to transfer the left atom 
into the first excited state of the right well without affect- 
ing the other one. This can be accomplished by changing 
the parameters Ui^ 2 (t) according to the following steps, 
which are illustrated in Fig. [5] 

(i) Between the times t — and t = t\ we raise very 
rapidly the minimum of the left well, such that its 
ground state crosses in energy the right well's first 
excited state. 

(ii) In the time interval \t\ , tg] we lower the central bar- 
rier down to a point where the atom can tunnel 
from left to right while at the same time we start 
to lower the left well. 

(hi) In the time interval [£2,^3] the barrier is raised up 
again while we continue to lower back the left well. 

(iv) During the time interval [ts , £4] we restore the initial 
potential shape. 

By doing steps (ii) and (iii) adiabatically the atom stays 
in the second excited state of the double well potential 
which is at t\ the ground state of the left well and at £3 
the first excited state of the right well. Thus the atom is 
transported from left to right. 

An effective transport procedure requires appropriate 
"pulse functions" U\(€) and u 2 (t) while the direction of 
the transport is governed by the parameters a, I. Let us 
assume for simplicity that the marker atom is located 
initially at a site kxj — (2j + l/2)ir, j 6 Z. Then the 
process shown in Fig. requires a — 1 and I = while 
for a further movement to the motional ground state on 
the right we have to set a = — 1 and I = 1 and to perform 
the same pulse functions backwards in time. For moving 
an atom from the ground state to the first excited state 
of the left neighboring well we have to set a — — 1 and 
I = 1 and to apply the forward pulse functions. In this 
case a further movement to the left ground state requires 
a = 1 and I = and the backward pulse function. Note 
that the abrupt changes of the phase (j> take place while 
U 2 is zero, i.e. when the corresponding lasers are com- 
pletely blocked off from the atoms. Every transport of 
an atom across the lattice can be divided into these four 
elementary processes. Since the pulse sequence is in all 
cases the same (except for time reversal) we can focus in 
the following on the example shown in Fig. [S] 

The feasibility of our quantum computing scheme de- 
pends on the time scale on which quantum operations can 
be performed. Clearly the latter is directly connected to 
the speed of the transport process. In this respect steps 
(i) and (iv) can be performed over much shorter times 
than step (ii) and (iii), which are limited for example 



by the energy difference to the other motional states. 
In order to examine adiabatic transport during step (ii) 
and (iii) it is thus necessary to study the instantaneous 
eigenenergies of an atom during the transport process in 
dependence on the parameters we can control, i.e. the 
pulse functions u\ t2 (t). 

Since U(x, t) is periodic, we can calculate the in- 
stantaneous eigenenergies and eigenfunctions of the sin- 
gle particle Hamiltonian by introducing Bloch functions 
ip^lx) — e lkx u^\x) with Bloch vector k and band in- 
dex n. In Fourier space the stationary Schrodinger equa- 
tion takes the form 



(q - fc) a u<"> (g, t) + J2 U(q - q', t)u[ n) (q, t) 



(15) 



where 



u^\q,t) = - P dxj« x u£\x,t), 



a-i Jo 



(x, t) = V e- iqx u k (q, t), q=—,neZ (16) 

* — * n ■ 



iik 



and 



U{q,t) = U a {t)S qfi +^^(8 2K , q + S 2K ^ q ) 

(17) 

where aj with i = 0, 1 are the lattice constants, i.e. 
ao = tt/k if U 2 — and 01 = 2it/k if U 2 7^ 0. By as- 
suming periodic boundary conditions the Bloch vector k 
gets quantized, i.e. k = 2nir/Mao, n = —M/2, . . . , M/2 
where M is the number of lattice sites. Given the func- 
tions Uk{x) we can furthermore construct Wannier func- 
tions which are localized at lattice sites xi, 



w 



(«) 



1 



'M 



E 



jk(x-Xi) 



<•""" l 'ui''(x). 



4°' 



(18) 



These functions are needed for example as initial states 
for solving the time dependent Schrodinger equation ifHJl. 
Since in our considerations the lowest bands are always 
practically flat, i.e. the Bloch states for a given band 
are approximately degenerate, the Wannier functions are 
also in good approximation eigenstates of the Hamilto- 
nian and thus dispersion of the wave packet during the 
time evolution is negligible. 

Equation i|15|) is a linear system of equations which can 
be solved numerically after truncating q at sufficiently 

high values. The instantaneous eigenenergies Ej^'(t) 
gained in this way are very important to find an adia- 
batic passage for the atom transport. An example for 
the band structure is shown in Fig. As can be seen 
from this figure, V is sufficiently large to be in a tight 
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FIG. 6: Band structures for two different values of u\,U2 and 
V = 100EV, cr = 1, I = 0. (a) Initial values of the lattice 
parameters (iii(0) = 112(0) = 0). (b) Band structure during 
the transport process (in = 0.99, ui = 0.13). The two figures 
correspond to the situations shown in Figs. [HJi and|KJ;. 



binding regime, i.e. the lower bands are flat and there is 
no tunnelling between different wells. 

For an efficient adiabatic transport the pulse functions 
have now to be chosen such that the corresponding en- 
ergies E^\t) behave in an appropriate way during the 
transport steps (ii) and (iii), i.e. one should for example 
avoid level crossings of the initial energy with the energies 
of other states. An example is shown in Fig. [7^. 

As already mentioned the Bloch states in the lowest 
bands are almost degenerate so we can restrict ourselves 
to a an arbitrary value of fc, e.g. k = 0. The upper solid 
line in this figure corresponds to the "path" of the atom 
to be transported initially located in the left well while 
the lower solid line indicates the path of the atom located 
in the right well. The initial depth of the potential wells 
is V = 100 E r leading to a "trap frequency" in a single 
well of about 20 oj r (u> r = E r /h is the recoil frequency). 
Keeping the height of the barrier constant we raise the 
minimum of the left well such that the ground state en- 
ergy of the left well crosses the first excited state energy 
of the right well [step (i)]. Then we proceed according to 
steps (ii) and (iii). By reducing the height of the central 
barriers the trap frequency decreases and by adjusting 
appropriately the pulse functions wc avoid that the solid 
lines cross the dashed lines. After time t% the original po- 
tential shape is restored [step (iv)]. The pulse functions 
Mi. 2^) used for this example and the corresponding time 
dependent physical relevant parameters J7i >2 (t), (p(t) are 
shown in Fig. [7jp and Figs. UJjd, respectively. 

The fidelities of the processes, given by 



F A = 



dx4> A (x,T)i>i a (x) 



(19) 



are numerically calculated by solving the time dependent 
single particle Schrodinger equation in position represen- 
tation i|14fl by using the Crank-Nicholson scheme [2{| 
where as initial state ip^ and final state i/>g we choose 
Wannier functions (|18|l which are located in the corre- 
sponding wells. The superscript A G {M, R} indicates 
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FIG. 7: Adiabatic atom transfer between lattice sites: (a) In- 
stantaneous eigenenergies iv,™ (t). The upper solid line cor- 
responds to the motional energy of the atom to be transported 
and the lower solid line corresponds to the register atom ini- 
tially in the right well (cf. Fig.^J. (b) Corresponding control 
parameters 111,2- (c) Amplitudes of the lattice components, 
(d) Phase difference between lattice components. The time 
axis is given in units of the recoil frequency u> r = E r /h. The 
time interval [ti,is] is not of the same scale as [0,ti], [t3,ti\. 



again the wavefunction of the atom to be transported 
(marker atom) and the atom which is supposed to stay 
located at its well (register atom). 

In case of the example of Fig. [7] we get F M = 99.91% 
for propagating the marker atom wavefunction and F R = 
99.98% for propagating the register atom wavefunction 
from t\ to £3 in a time T — £3 — t\ = 2Q/oj r . In the case 
of Rubidium (uj r = 2ir x 3.8 kHz) this would correspond 
to a time T = 0.8 ms and for sodium (ui r = 2ir x 25 kHz) 
we would have T = 130 /its. For the optical superlattice 
described in [24| a laser power of P = 3 mW was suffi- 
cient to create a maximal potential depth of 217; < 60-Er 
for 87 Rb. Keeping the ratio of laser intensity and detun- 
ing constant we have £7; ~ VP. For a potential depth of 
roughly Ui < 50Er which is required in the above exam- 
ple (see Fig. Et) we can estimate the required maximal 
laser power to be merely P«8-9 mW. 
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C. Pulse optimization 

If we relax the constraint of adiabatic transport dur- 
ing step (ii) and (iii) the process described in the previous 
subsection can be significantly accelerated. In this case 
the pulse sequences have to be engineered in a certain 
way which can be done by using quantum optimal con- 
trol techniques as detailed, e.g., in [2^, [53, [2^1 . Thereby 
the evolution of a quantum system governed by a set of 
control parameters (in our case these are the functions 
ui, 2 (t)) is tailored to reach a pre-determined target state 
tpRn with optimized fidelity within a specific time T. For 
notational convenience we will denote in the following the 
time t = ii as t = and t = t 3 as t = T. 

The basic idea is to minimize the infidelity of the pro- 
cess with the constraint that the Schrodinger equation 
has to be fulfilled. This amounts to find the station- 
ary point of a functional, leading to a set of equations for 
the wavefunction and auxiliary states \ A which are intro- 
duced as Lagrange multipliers. In our case this functional 
takes the form 

C^ M ^ R A M A R ,X M ,X R ,u u u 2 ) = 



E 

Ae{MM} 



dx^£ n (x)il> A (x,T) 



2Ro| J dt J dxx A (x,t)(ip A (x,t) 



(20) 



with 



h 2 d 2 

H{ Ul (t),u 2 (t)) = - — — + U(x,u 1 (t),u 2 (t)), (21) 
2m ax 

where the last term in Eq. (|21|l is the potential © , which 
depends on the control parameters. As can be seen from 
Eq. I|20l) we are looking for a minimum of the sum of the 
infidelities of the process for the marker atom and the 
register atom. Setting the derivatives with respect to the 
arguments of C equal to zero leads to the following set of 
equations, 



ihip A (x, t) = H( Ul (t),u 2 (t))i> A {x, t), A e {R, M}, 



(22) 

ihx A (x,t) =H( Ul (t),u 2 (t)) X A (x,t), A&{R,M} 

(23) 

with conditions 

iP A (x, 0)=^ A dx), (24) 
X A (x,T) = ^ A n (x) J dx^£ n (x)^ A *(x,T), (25) 

and 

= -2Im K A {u x {t),u 2 {t)), j = 1,2 (26) 

Ae{M,R} 



with 
Kf(ui(t),u 2 (t)) 



dxip A (x,t) 



dll{x,ui(t),u 2 (t)) A * 



hduj (t) 



X A (x,t). (27) 



These equations are the basis of the optimal control al- 
gorithm which minimizes C with respect to uj . Thereby 
we solve Eq. 1(22(1 and 1)23(1 numerically by introducing a 
discretized time axis with time step At and by using the 
Crank-Nicholson scheme. For the sake of completeness 
we briefly describe the algorithm we use here. The fol- 
lowing procedure, called immediate feedback control, is 
guaranteed to give a fidelity improvement at each itera- 
tion [23. 

The Schrodinger equations 122(1 are integrated from 
t = to t = T leading to ip A (x, T) with an initial guess 
for the control parameters uf\(t). At this point an iter- 
ative algorithm starts during which the controls Ui 2 (t) 
are updated. 

Let us assume that we are in the n-th iteration. Tak- 
ing the controls u^, ' 2 (t), Eqs. have to be solved back- 
wards in time, i.e. from t = T to t = 0, with "end values" 
(|25(l which can be interpreted as the part of tl) A {x, T) that 
has reached the objective. Given the solutions x A { x ,fy 
the functions x" 4 (x,t) and ip (x,t) (with initial condi- 
tions 1(24(1 ) are now again evolved forward in time while 
the control parameters are updated during each time step 
according to 



(n+l) 



(t) 



(*) 



\(ty 



-Im 



£ Kf{u^{tU 2 n \t)) 



Ae{M,R} 



(28) 

During the forward evolution x(cc,i + At) is calculated 
using the controls u X2 (t) while ip(x,t + At) is evolved 

according to uj" 2 +1 ^ (t). The weight X(t) is used to enforce 
fixed initial and final conditions on the control pulses. 
Given these solutions we go on with the next iteration. 

The results of such a calculation after 135 iterations 
are shown in Fig. [S] which corresponds to the situation of 
Fig- El Fig.|SjL shows the control parameters for transfer- 
ring the marker atom to the first excited state of its right 
neighboring site while keeping the register atom at its ini- 
tial location. Fig. [SJdc shows the corresponding physical 
parameters. As starting values uf\{t) we took the pulses 
of the adiabatic example (see Fig. |7J). The use of the 
optimized pulses leads to a reduction of the transport 
time down to T3 = hh/En, which corresponds to 200 /is 
for rubidium and to 32 /is for sodium with a fidelity of 

F M = pR = 99.99%. 



IV. COHERENT RESONANT COLLISIONS IN 
A TRAP 

The coupling scheme we are proposing can be imple- 
mented either in dipole-force potentials, like optical lat- 
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(ii) 


(iii) 


(iv) 
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FIG. 9: Interaction of two atoms in an external trap. The two 
curves schematically describe the Born-Oppenheimer poten- 
tial for two scattering channels around a Feshbach resonance. 



FIG. 8: Non-adiabatic atom transfer between lattice sites: (a) 
Control parameters ui,2- (b) Corresponding amplitudes of 
the lattice components, (c) Phase difference between lattice 
components. The time axis is given in units of the recoil 
frequency ui r = E r /h. The time interval [ti,ts] is not of the 
same scale as [0,ii], [t3,i4]. 



tices l3(M . or in static electromagnetic traps, like atom 
chips 21] . Performing gate operations as described in 
Sec. IIVI requires a strong molecular interaction between 
register atom and marker atom. Atoms can be coupled to 
molecular states either by means of Feshbach resonances 
[22^ or through Raman photoassociation laser pulses [3l| . 
For the sake of concreteness, we focus here on Feshbach 
resonances in optical lattices - however, all of our argu- 
ments can be adapted, e.g., to Raman photoassociation 
on atom chips. We consider 87 Rb atoms trapped in a 
two-component optical lattice (see, e.g., [24j). 



A. Feshbach resonances in confined geometry 

A schematic picture of the Born-Oppenheimer poten- 
tial describing their interaction in the relative coordinate 
r is shown in Fig. |5J Negative values of the excitation 
number v label bound molecular eigenstates of the dimer 
system, while positive v values denote unbound trapped 
two-atom states. Such a potential exists for each col- 
lision channel \(3), corresponding to the relative-motion 
and hyperfme angular momentum quantum numbers of 
the two colliding atoms. Feshbach resonances occur when 
a bound state \np) (v = —2 in the example shown) crosses 
the dissociation threshold for a state having the same 
quantum numbers | 22| while changing an external mag- 
netic field B. Close to resonance, the scattering length 



varies as 

W = i6j ( 1_ ^)' 

where Ab g is a non resonant background scattering 
length, Bp is the resonant magnetic field, and is the 
width of the resonance. The resonance energy varies al- 
most linearly with the field 

e p (B) = sp(B - Bp), (30) 

with a slope sp. We are interested in the dynamics of 
such a system in a confined geometry. Following |32l | , we 
shall model it by the effective Hamiltonian 

Hp = sp{B)\np){np\+Y J («M«) (« \ + V v\ «> (np\+h.c), 

V 

(31) 

where the |f)'s are the trapped relative- motion atomic 
eigenstates of an isotropic harmonic oscillator trap hav- 
ing frequency v. The couplings to the resonance are 

Kf = 2hvy/\/4^+3 a bg 6p/ir (32) 

with abg = Ab g y/mv/h, 8p = Apsp/(hu). In a differ- 
ent geometry, for instance in an elongated trap charac- 
terized by a ratio 7 between the ground level spacings 
in the transverse and in the longitudinal potential, the 
couplings can be calculated by projection on the corre- 
sponding eigenstates (see Appendix lA"|l . Accurate values 
for the resonance parameters Ap and Bp, as well as for 
A b g, are now available from both theoretical calculations 
and recent measurements j3^|. 

The possibility of controlling the resonance energy via 
an external magnetic field, as described by Eq. I|3U|) , pro- 
vides a straightforward way to steer the interaction be- 
tween the atoms. Not only can the scattering length be 
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FIG. 10: Adiabatic (solid) and diabatic (dashed) energies of a 
6-level expansion (molecular resonance level plus 5 trap levels 
of an isotropic 100 kHz trap) with a linear ramp of magnetic 
field with slope 5 G/ms. The adiabatic energy levels show a 
set of avoided crossings. We use the Feshbach resonance level 
near 100 mT. The coupling matrix element Vbis 0.884 hv for 
this trap. 



B. Choosing qubit logical states 

We will identify our qubit logical states with the clock- 
transition states 

|0) = |F = 1,tof=0), |1) = \F = 2,m F = 0) (33) 

and the auxiliary state as 

\x) = \F = l,m F = l). (34) 

The main advantage of this choice is that the qubit states 
are not sensitive to the magnetic field, and hence not 
subject to decoherence due to its fluctuations. 

The level scheme for a single atom is shown in Fig. 1111 
When we consider two atoms, the relevant level scheme is 



-1 



|1> 



+1 



+ 2 



F = 2 



10} 



\x) 



F = l 



varied over a significant range - the atoms can also be 
adiabatically coupled into a molecular state. An example 
of this sort of process is shown in Fig.^| Here, the eigen- 
values of the interacting Hamiltonian are shown for a six- 
level model including the five lowest unbound trap states 
plus the resonant state. The latter is ramped across 
threshold by applying an external magnetic field having a 
linear dependence on time. Both the so-called "diabatic" 
energies (i.e. those obtained by neglecting the couplings 
to the resonance) and the adiabatic ones (i.e. the actual 
eigenvalues of the full coupled Hamiltonian) are plotted 
against time for a certain ramp rate. The important 
point to notice is that the ground state of the relative 
motion is adiabatically connected to the resonant state. 
Therefore, if the atoms are prepared in their relative- 
motion ground state and the resonance state is ramped 
across threshold from above, the atoms are transferred 
into the bound state, whose energy depends on the mag- 
netic field - and the process is actually reversible. This 
mechanism has been used for the creation of molecules 
in ultracold gases |H |H EHH |H HJ , and is very 
relevant in the present context of quantum information 
processing due to its inherent state-dependent nature. 
Indeed, the coupling to a specific resonant state |ng) 
is only effective for a particular entrance channel 
while in general all other combinations of atomic hyper- 
fine states (that is, of logical qubit states in our case) 
will be unaffected by the resonance. Thus the resonance- 
induced energy shift will cause a two-particle phase to 
appear only for that particular two-qubit computational 
basis state. We will see in the next Section how to use 
this effect in order to achieve a desired C-phase gate. 



FIG. 11: Internal level scheme for a single atom: qubit and 
auxiliary states. 

described by Fig. 1121 if they occupy their relative-motion 
ground state. Appendix iBl shows that this is indeed the 
case for two bosons stored in the ground and first axial 
excited state of a cigar-shaped harmonic trap. In Fig,IT2l 




FIG. 12: Internal level scheme for two coupled atoms: levels 
involved in single- (left) and in two-qubit operations (right). 
The Raman transitions fh (fb) used for single- (two-) qubit 
operations are shown. 

the resonant levels \np) are also shown, which are used to 
induce energy shifts for the purpose of gate operation as 
described in the next Section. Indeed, in a confined ge- 
ometry, the coupling to such molecular states can induce 
dressing of the trapped eigenstates with a half-splitting 
Ae' 3 - controllable, by varying the external field, up to 
a maximum value equal on resonance to the interaction 
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FIG. 13: Dependence of the scattering length of Rb on the 
external magnetic field for collisions in channels 1 00) and |Cte). 



strength V@ -, as shown in Fig. El f° r the collisional 
channels \j3) = |00), ]0ar>. 

Our calculation with a realistic molecular interaction 
potential yields, among others, two resonances at Bq x = 
386 G and B ao = 407 G, having widths A 0x = 5.7 mG 
and Aoo = 16 mG, as shown in Fig. ^5] These will be 
employed in the following to effect logical gate operations. 



cific sites where the single-qubit operation takes place 
can be selected even if the addressing laser cannot re- 
solve them spatially from neighboring sites. Moreover, 
the two-atoms state remains always factorized, whence 
possible magnetic field fluctuations, affecting the state 
|x) (unlike |0) and |1)) will yield only a global phase. In 
a rotating frame the Hamiltonian which describes this 
system takes the form 

H =6\lx)(lx\ + e 0x (B)\M)(M\ 

+V 0x (\0x)(M\ + \M)(0x\) + ^(|0x)<lx| + \lx)(0x\) 

(35) 

where <5 is the Raman detuning of two co-propagating 
Raman lasers, Vq x is the coupling between the molecu- 
lar state \M) and the dissociated state |0x) and so x (B) 
is the energy of the molecular state (we set the energy 
of the | Ox) state to zero). At the beginning of the op- 
eration the lasers are switched off (i.e. f2i = 0) and the 
external magnetic field is adiabatically tuned to the Fes- 
hbach resonance, i.e. sq x {B) — > So x (B 0x ) = 0. This leads 
to a splitting of the two particle state |0x) in two new 
eigenstates, 



\sf) 



1 

71 



(|0x)±|M)). 



(36) 



These states are indicated in the left hand side of Fig.HTl 
(here we have Ae 1 



Ox 



yOx 



since we are on resonance 



For finite laser power, in this basis the Hamiltonian takes 
the form 



V. QUANTUM OPERATIONS 

Let us now examine in detail how to use the features 
described above in order to perform quantum computa- 
tion in our system. Preparing all atoms in an initial state 
by optical pumping requires no single-qubit addressing 
and can be performed with standard techniques. As the 
next step, according with the above discussion, perform- 
ing a specific algorithm will require a certain pattern of 
marker atoms. These can be prepared either in a periodic 
fashion, by means of a supcrlattice tuned to the appro- 
priate transition, or in an ad hoc lattice region, spatially 
separated from the one where computation has to take 
place, to be subsequently loaded into the latter via the 
transport mechanism detailed above. 



A. Single-qubit gates 

In the single-qubit case, the relevant resonance field is 
Bq x , while the lasers couple |lx) with the lower dressed 
state connected to |0x) with an effective Rabi frequency 
Sli as in Fig. ^Jb. The process is resonant only if the 
marker atom in state |x) is present. In this way, spe- 



H =tf|lx)(lx| + V£ x (\Sf ){Sf\ - |S te )(Sr|) 



2V2 



(|^) + |5 te ))(lx| + |lx)((^| + (5 to |) 

(37) 



The Raman detuning is set to 5 = — Vq x and if 
£1 /Vq x « 1 we can project out the state \S^). This 
yields the two level Hamiltonian 



H 



y te (|ix)(ix| + |^)(^'|) 

^(|S te )(lx| + |lx)(S to |) (38) 



i.e. if we finally tune the magnetic field out of the Fesh- 
bach resonance again we get the transformation 

|0x) ->-cos(fii/2\/2£)|lx) -isin(fii/2V2i)e i¥> |0x), (39) 

|lx) ->-isin(Q 1 /2\/2t)|lx) + cos(Q 1 /2v / 2t)e i¥ '|0x). 

(40) 

In this expression we included a phase ip which is the (ad- 
justable) phase accumulated during the adiabatic ramp- 
ing process of the magnetic field. 
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B. Two-qubit gates 

In the two-qubit case, we take the marker atom to be 
in a state of the logical subspace spanned by |0) and |1). 
This time, the field is ramped across -Boo i an d the Raman 
lasers couple for a time r - with Rabi frequency il 2 ~ the 
lower dressed state to the degenerate two-atom levels 1 01) 
and 1 10) (Fig. 1121c). In a rotating frame the Hamiltonian 
can be written as 

H =tf(|01)(01| + |10)(10|) + 25|11)(11| 



e 0Q (B)\M)(M\ 
2 



■(|00)<01| 



(|M)(00| + |00)(M|) 
|00)(10| + |11)<01| + |11)<10| 



I- h.c. 
(41) 



The notations are the same as in Eq. J3HJ- We perform 
now the same procedure as in the case of the single qubit 
rotation, i.e. we tune adiabatically the magnetic field to 
the Feshbach resonance, i.e. eoo(-Boo) = while Q2 = 0. 
The Hamiltonian with diagonalized molecular part reads 

H =<5(|01)(01| + 1 10) <10|) + 2<S|11)(11| 

+ v 00 (\s + Q )(s° + Q \-\s 00 )(s 00 \) 

Mr 



2V2 
|11)(01 



\s 



(III 



> + |S°°))<01| 
11)(10| + h.c. 



(|fl?> + |5?»(10| 



where 



\S° ± °) 



V2 



(|00)±|M)), 



(42) 



(43) 



These states are shown on the right hand side of Fig. ITT1 
(now we have Ae 00 = Vq°). Taking the Raman detuning 
to be S = -V$° amounts to the fact that (if 9, 2 /V§° < x ) 
the states |11) and jS 1 ? ) are effectively decoupled from 
the remaining three states. Projecting out the uncoupled 
states the effective Hamiltonian for the remaining three 
level system then takes the form 



H 



V 00 (\S m )(S m \ 
Q ° (|S 00 ) (01 1 



2V2 



|01)(01| + |10)(10|) 
|5 oo )(10| + |01)(5 oo | 



|10)(S 00 |). 
(44) 



If we introduce the vector notation \ip) <-> 

((S ,Q _°|?/'), (01 (10| V}) an d disregard global phases 
the time evolution operator of this system can be written 
as 



, / 2c(t) -is(t) -ia(t) 
U(t) =- -is{t) c(t) + 1 c(t) - 1 
1 \-is(t) c(t)-l c{t) + l J 



(45) 



with c(t) — cos(f22^/2) an d s(t) — cos(51 2 </2). If we 
apply a Raman pulse of duration r = 2(2n + l)7r/S^2 



PI 
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FIG. 14: Two-qubit (C-phase) gate operation: resonance en- 
ergy £00 (top; off-resonant before and after gate operation), 
overlap between evolved and initial state (center), accumu- 
lated phase (bottom). 



and finally tune the magnetic field out of the Feshbach 
resonance again, we get the following truth table for the 
operation: 



|00) - 


-> -e^|00), 


|01) - 




|10) - 


--|io>, 


111) - 





(46) 



where we included again the phase tp now accumulated 
by state 1 00} during the ramping process due to the 
interaction energy shift, whose value can be adjusted 
by controlling the magnetic field. For ip — 2t: and 
Vo°°t = (2m + l)7r - which imposes a commensurability 
condition ^2/^0°° = 2(2n+l)/(2m+l) between the Rabi 
frequency and the Feshbach energy shift -, a swap oper- 
ation is performed. Besides being an essential ingredient 
for entangling gates between distant atoms as detailed 
in Sec. ^ such a swap operation can greatly help in the 
task of non-destructive qubit readout. To this aim, the 
quantum state of an atom to be read out at the end of 
a computation could be simply swapped onto a marker 
atom to be subsequently transported to a different lattice 
region where measurement can take place without phys- 
ically disturbing the register atoms, which can be later 
re- used for logical operations. 

On the other hand, if no Raman lasers are present 
and ip — 7r, a C-phase gate between register and marker 
atom is obtained. A two qubit gate between distant reg- 
ister atoms can be realized as described in Sec. |H] Note 
that laser addressing of single qubits is never required 
throughout the procedure. 

The magnetic ramping process can be even performed 
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non-adiabatically, provided that all population is finally 
returned to the trapped atomic ground state. This can be 
accomplished via a quantum optimal control technique in 
analogy with the above discussion for the transport pro- 
cess. The control parameter in this case is the resonance 
energy eoo, which can be adjusted by varying the exter- 
nal magnetic field. Care has to be taken in optimizing 
not only the absolute value of the overlap of the final 
state onto the goal state, but also its phase <p. Fig. ^] 
shows the optimization results for a 100 kHz trap with 
a ratio of u±/v = 10 between the trap frequencies v_\_ in 
y, z-direction and v in x-direction. The final infidelity is 
about 2 x 10 -5 in this case. 



VI. CONCLUSIONS 

When it comes to using neutral atoms for the pur- 
pose of quantum information processing, besides the well- 
known general criteria formulated by D. DiVincenzo |4l| . 
the fulfillment of various practical requirements, specific 
to atomic implementations, can make a difference on the 
road to experimental realization. For example laser ad- 
dressing of single qubits, though being theoretically triv- 
ial, is limited by diffraction, imposing a lower bound on 
the actual spacing between qubits. Furthermore perform- 
ing gate operations in state-dependent potentials cre- 
ates entanglement between internal and external degrees 
of freedom, which in turn is prone to decoherence, as 
random fields typically affect differently the two logical 
states. The same is true for internal-state entanglement, 
if the qubit states are chosen with different Lande factor 
and, unless the latter vanishes for both states, they will 
be sensitive to magnetic-field fluctuations. 

In this paper, we introduced the concept of "mark- 
ing" qubits via molecular interactions which allows for 
relaxing a number of these constraints for neutral-atom 
quantum computing. We have presented a scheme that 
enables quantum gates and information transport in a 
quantum register, even though requiring neither single- 
site addressing by externally applied fields nor state- 
dependent external potential. Moreover, qubit states 
with the same (even vanishing) Lande factor can be em- 
ployed; and the overall speed can be of the order of the 
inverse atomic trapping frequency. We have shown how 
this scheme can be implemented in two-component opti- 
cal lattices, whereby the mechanism used to mark atoms 
is the molecular interaction responsible for Feshbach res- 
onances, which are currently a subject of intense experi- 
mental research in the field of cold atoms, where molecule 
formation via control of Feshbach resonances has been re- 
cently achieved IH S E3 S |M E3 • In other words, 
our proposal relies on techniques that are presently being 
developed, and represents therefore a feasible candidate 
for the implementation of quantum information process- 
ing with neutral atoms in optical lattices. 

Finally, the analysis presented here is limited to one- 
dimensional systems, basically with a single marker 



atom. Further conceptual development is possible, for in- 
stance in exploring the interplay between several marker 
atoms on the same lattice, or the extended flexibility 
given for instance by higher-dimensional geometries; this 
will be the subject of future investigations. 
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APPENDIX A: DYNAMICS IN A 
CIGAR-SHAPED TRAP 

The normalized eigenfunctions of a 3D harmonic oscil- 
lator in spherical coordinates are 



*«AnA r '">W y nlT{£ + 3/2) 2 



r e 



!F!(-n,£ + 3/2, aV)y/™* (<?,</>) (Al) 



where a — y/mv/h, T(n) is the Euler Gamma func- 
tion, \Fi(a, 6, z) is the Kummer confluent hypergeomet- 
ric function, and Y™ e (6, </>) are the spherical harmonics. 
The normalized eigenfunctions in cylindrical coordinates 
(we assume the same frequency v in the longitudinal di- 
rection, and a transverse frequency v±_ a factor 7 higher) 
are: 

- vm7m> — x ,A2) 

H n {az) xFti-e^me] + 1,70V)/ 



A (-iM~i) 3 r(K| + » + j + i) 



where H n {x) are the Hermite polynomials, and (a)j = 
r(a + i)/T(a) is the Pochhammer symbol. We are inter- 
ested in s-wave scattering processes, so we restrict our 
analysis to the eigenstates with i = mi = and obtain 



cyl 



IC P nn) =0; 



2r+l,0,0 I ~*w,0,Q/ ~ u ' (A3) 



i=0 j=0 

1F1 (j + \ 



(2z + l)!(2j + l)! 



2i'~J~2'-' T 2' 7+I 



(v- j)l(w-i)l(j + l)*+J"+3/2 
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The coupling matrix elements in a cigar-shaped trap 
with anisotropy factor 7^1 are computed as 

Vf (7) = E<^ool^ P oo)^(7 = 1), (A5) 



where the spherical matrix elements V w (^f = 1) are given 
by Eq. p?f. 



l*i,o) 



II), |0> 2 + 10), |1) 2 ^1)1^0)2 + ^0)1^1)2 



V2 



V2 



|1>! |0> 2 + |0>X |1> 2 1 , i 111 

— - — ^ — -m)xm)x, 
1^1)11^0)2 + 1^0)11^1)5 



1)1 |i) 2 



V2 

1>, |1> 2 |^l>^|^0>^- 



(B3) 



(B4) 



APPENDIX B: CONDITIONAL LEVEL SHIFT IN 
A QUASI-1D TRAP 



Let us consider the three-dimensional state \^s,m s ) 01 
two spin- 1/2 bosons in a harmonic trap. In the x direc- 
tion, one particle is in the trap ground state IV'o), and the 
other in the first excited state l^i)- The transverse state 
is the ground state |^j_) for both particles. X and x are 
the center-of-mass and relative coordinate. Denoting the 
i th particle's state by |-),-, the symmetrized states can be 
written as 



l*o,o) - 



|1>! 10)2-10), |1) 2 |^l)l|^0)2 -|^0)l|^l)2 



V2 

|1),|0) 2 -|0),|1) 2 



l*i,-i) = |0)x |0) 



V2 

- \ipo)x\ipi} x ; 
1^1)1^0)2 + 1^0)1^1)2 



V2 

1 1°) 2 \A)x\iPo) x ; 



(Bl) 



(B2) 



When we apply a static external magnetic field corre- 
sponding to the Feshbach resonance for the |00) channel, 
the interaction only affects the state |\l/i._i), dressing it 
with a splitting 2Vq 00 (7) that can easily be of the order of 
Hv assuming a ratio Vx/v = 10 of the trap frequancy v±_ 
in y, z-direction and v in x-direction. This means that 
the state |0), |0) 2 can be discriminated spectroscopically, 
allowing for different kinds of gate operation as described 
in the text. 
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